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' Abstract 

The capacity region of a channel consists of all achievable rate vectors. Picking a particular point 
in the capacity region is synonymous with rate allocation. The issue of fairness in rate allocation is 
^ addressed in this paper. We review several notions of fairness, including max-min fairness, proportional 

, fairness and Nash bargaining solution. Their efficiencies for general multiuser channels are discussed. 

We apply these ideas to the Gaussian multiple access channel (MAC) and the Gaussian broadcast 
channel (BC). We show that in the Gaussian MAC, max-min fairness and proportional fairness coincide. 
^ ' For both Gaussian MAC and BC, we devise efficient algorithms that locate the fair point in the capacity 

. region. Some elementary properties of fair rate allocations are proved. 

1 Introduction 

, There are several fairness criteria in rate allocation. The simplest one is to mandate that all users have 

I the same data rate and maximize this common data rate within the capacity region. This allocation is 

O ' equitable and maximize the data rate of the worst user. However, it often does not effectively utilize 

system resources. 

Another criterion is called max-min fairness. It relaxes equity and allows increasing the rates of some 
■ users without lowering the minimum data rate in the system. Under such an allocation policy, nobody can 

5^ I be benefited by worsening anybody who has lower data rate. It can be rephrased as follows. If we take 

some resources from a wealthy user to a relatively poor user (without reversing the order of wealthiness), 
the resulting allocation is considered fairer. Such operation is called a Robinhood operation. We say that 
an allocation is max-min fair if no Robinhood operation is possible without violating feasibility. 

The max-min fairness is considered quite stringent. Any decrease of rate of a user with low data 
rate cannot be compensated by increasing the rate of any user with higher rate, no matter how small the 
decrease is. Kelly considers logarithmic utility function and proposes proportional fairness for network flow 
control [1]. Roughly speaking, we say that a rate allocation is proportionally fair if any adjustment will 
decrease the sum of percentage change over all users. A framework of optimization is thereby introduced 
using Lagrangian technique. Both max-min and proportional fairness are popular criteria for flow control 
in unicast and multicast networks [SJ |3] . 

Another classical notion of fairness is the Nash bargaining solution coined by Nash in the 50's 4J. In 
fact, proportional fairness is its special case. The Nash bargaining solution is a standard tool in cooperative 
game theory, and is applied widely in network resource allocation. For example, see [5] for application to 
orthogonal frequency division multiple- access networks. 
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Yet there is another fairness criterion based on the theory of majorization. We wiU give formal defini- 
tions of all the above criteria in the next section. 

In this paper, the problem of picking a point in the capacity region of Gaussian MAC and BC according 
to some fairness criteria is considered. Some works have been done on Gaussian MAC. The sum rate at 
the max-min fair allocation is characterized in '6^. The polymatroid structure of the capacity region 
of Gaussian MAC is exploited in [71 EI, and algorithms for locating the max-min fair point are found. 
Implementation issue is addressed in j9]. 

The organization of this paper is as follows. In Section [21 we give precise definitions of fairness, and 
review the theory of majorization and Schur- convexity. The results are summarized in Section [3l and the 
details for the Gaussian MAC and BC are in Section [H and [5l respectively. The appendix contains some 
proofs of the theorems in Section [H 



2 Fairness, Majorization and Schur- convexity 

The set of all users is denoted hy Q = {1, . . . , K}. 

In this section, we review several fairness criteria, and the theory of majorization. We will use the 
symbol TZ to represent capacity region, which is assumed to be a closed and convex set throughput this 
section. 



An allocation is symmetric if every user has the same data rate. Symmetric capacity is the maximal 
sum rate of all symmetric allocations, 

CsymCR-) := ma,x{Kr : {r,r, . . . ,r) e TZ}. 

An allocation is called max-min fair if we cannot increase the rate of user i without decreasing rj 
for some rj < ri, while maintaining feasibility. At the max-min fair allocation, no user can increase the 
data rate without compromising users with lower data rate. Formally speaking, a rate allocation r-^^^ is 
max-min fair in TZ if for any r G TZ such that rf^^^ < ri for some i, then we can find j € ^ such that 
rj < rj*^*^ < rf'^^'^ . The sum of rate at the max-min fair allocation for capacity region TZ is called the 
max-min capacity and is denoted by Cmm{TZ). 

Proportional fair (PF) rate allocation {rf^)i=i,,,K is the data rate allocation that maximizes 

K 

1=1 

The proportional fair capacity is the corresponding sum rate, 



Since the capacity region TZ is closed and convex, and the log function is concave, the maximization is 
well-defined. Another characterization of the proportional fair allocation is 



K 



i=l 



PF 



for all point r in the capacity region TZ. 

In the Nash bargaining solution, there is a notion of disagreement point, which is the default operating 
point if the users fail to reach any agreement. User i will not accept any data rate lower than di. The rate 
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Figure 1: The point d is the disagreement point. The shaded area is the acceptable rate allocation. 

allocated to user i in the Nash bargaining solution should be larger than or equal to di (Fig.[T]). If we are 
given a disagreement point d in the capacity region, the Nash bargaining solution maximizes 

K 

^log(r, - di) 

1=1 

over all points in the region 

{r e 7^ : r,> di Vi}. 

The Nash bargaining solution satisfies several desirable properties. See fTU| for details. It is obviously 
identical to the proportional fair solution when the origin is chosen as the disagreement point. 
For a vector x = {xi, . . . xk) G we denote the components in nondecrcasing order by 

X[l] < X[2] < ■■■ < X[K]- 

We say that vector x is majorized by vector y, written as x ^ y, if for k — 1, . . . , K — 1, 

k k 

i=l i=l 
K K 

i=l i=l 

In other words, if we sort the components of x and y in nondecrcasing order, the cumulative sum of the 
components of x is larger than the corresponding cumulative sum of y. 

Majorization induces a partial order that measures dispersion. It is known as the Lorenz order in 
economics, and is used for comparing income distributions. When x ^ y, we say that the distribution 
according to x is less spread out, and is thus fairer than that of y. A canonical example is that the vector 

is majorized by any vector in E;^ whose components sum to 1. 
A function / : M.^ ^ R is called Schur convex if 

/(x) < /(y), whenever x ^ y. 

If the inequality above is reversed, then we say that the function / is Schur concave. A class of Schur-convex 
functions is constructed using the following lemma. 
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Lemma 1. //0 : M ^ M is convex (concave), then the function 

K 

i=l 

is Schur- convex (Schur- concave). 

A useful criterion for Schur-convexity is as follows. 

Lemma 2 (Schur's criterion). Suppose that F : M.^ — > K is differentiable and symmetric, meaning that 
F{xi, . . . , Xk) = ^(2^77(1) ; ■ • ■ J X-rrlK)) for any x and permutation tt 0/ {1, ... , K}. Then 

1. F is Schur- convex if {xi — Xj) ^^^(x) — ^^(x)^ > for all i and j , 

2. F is Schur-concave if {xi — Xj) ^^(x) — ^^(x)^ < for all i and j. 
See for more on the theory of majorization and Schur-convexity. 

3 Capacity-Fairness Tradeoff 

Usually there is a tradeoff between sum rate and fairness. The next theorem illustrates such a tradeoff. 
It shows that the requirement of symmetric fairness is more stringent than proportional fairness in the 
sense that the symmetric capacity is always less than or equal to the proportional-fair capacity. We use 
the notation Csum for the sum capacity, defined as 

K 

Csum{'R-)= max 'S^Ji. 

2 — 1 

Theorem 3. For any convex region TZ, 

C syjn 

{n)<CMM{n)<Csum{'R). 

Equality holds in the first inequality in the first line only if the PF allocation is symmetric. 

Proof. All inequalities are obvious except the the first inequality in the first line, i.e., CsyryiiTV) < CppiTi-)- 
Let (ro, . . . , To) be the maximal symmetric rate allocation, and (rf ^, . . . , r^^) be the proportional fair 
rate allocation in region TZ. We have 

i=l i=l 

The first inequality comes from the defining property of proportional fairness, and the second is the 
AM-GM inequality. Thus, 

1=1 

When equality holds, then we must have rf^ = rl'^ = . . . — . This proves the first inequality. □ 

In general, the max-min capacity may or may not be larger than the proportional fair capacity. How- 
ever, in Gaussian MAC, both max-min and proportional fair capacity achieve the sum capacity. In fact, 
the max-min fair and proportional fair rate allocation in Gaussian MAC coincide. It is hence not necessary 



4 





Symmetric 


Max-min and PF 


Max. sum rate 


V 


[0,1] 


1 


1 


C» is Schur concave in P 


Theorem [S] 


trivial 


trivial 


Pi< ...<Pk ^rl< ...<r*j, 


trivial 


Theorem 1151 


N/A 


P ^ P ^ r* ^ f * 


N/A 


Corollary [T7] 


N/A 



Table 1: (Summary of results for Gaussian MAC) Capacities as a function of the power constraint vector 
P are Schur concave. For max-min fair or proportional fair, the rates are in increasing order if the power 
constraints are sorted in increasing order. If the power constraints are more spread out, so does the 
associated rate allocation. 
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PF 


Max. sum rate 


V 


[0,1] 


[1/^,1] 
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C* is Schur convex in N 


Theorem!^ 


Conjecture 1331 


trivial 


Ni < ... < Nk ^rl > ...>r*j^ 


trivial 


Corollary [30] 


trivial 


Ni<...<Nk^pI<...<P*k 


Theorem [25] 


Theorem [32] 


false 



Table 2: (Summary of results for Gaussian BC) Capacities as a function of the noise vector N are Schur 
convex. When the noise powers are sorted in increasing order, then the data rates are in decreasing order, 
and the power are in increasing order. 

to distinguish between max-min and proportional fairness in Gaussian MAC. Meanwhile, in Gaussian BC, 
max-min capacity is the same as symmetric capacity. 

Given a rate allocation r, we define its efficiency, i], as the ratio between the sum rate at r and the 
maximal sum capacity. Obviously, this is a number between zero and one. For proportional fairness, we 
have the following lower bound on rj. 

Theorem 4. Let TZ be a convex region in R;^, the efficiency for proportional fairness, rjpp, is lower 
bounded by l/K. 

Proof. By the definition of proportional fairness, the region TZ is contained in the polyhedron V defined 

by 

and ri > for all i. Therefore, Csum{T^) is not larger than Csum{P), which is equal to K maxi rf^. Hence, 

Csum (7?.) K maxi rf ^ K 

□ 

We will show that in Gaussian MAC, rjpp is exactly 1; it achieves the maximal value. In Gaussian 
BC, the lower bound l/K is attained, i.e., inf Cpp/Csum — l/K with the infimum taken over all K-usei 
Gaussian BC. 

The major results are summarized in Table [T] and [2] For both Gaussian MAC and BC, we also devise 
efficient algorithms that compute the fair solutions. 
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Figure 2: Capacity Region of a 3-user Gaussian MAC channel 

4 Multiple- Access channels 

In a scalar Gaussian MAC with K users, the received signal is 

K 
i=l 

where Xi is zero-mean Guassian with variance at most Pi and Z is Gaussian noise with power N. Let C 
denote the Shannon capacity formula 

C{x) = ^log{l + x). 
The capacity region of a scalar Gaussian MAC is [12] 

|r e : ^ i?^ < C(— P») , for all S <Z fl^ 

ies i&s 

The faces of the capacity region are hyper-planes in the form X^jes ~ ^'-"^ some subset S" of 17 and 
constant c. 

In vector Gaussian MAC, the received signal is 

K 



Y = ^ X,s, + Z. 

1=1 

The Si's are unit-norm column vectors of length L, Xi is a Gaussian random variable with zero mean and 
variance at most Pi, and Z is the Gaussian noise vector with zero mean and covariance matrix iVI, where 
I is the L X L identity matrix. The capacity region of a vector Gaussian MAC is [T3] 

{r e : r, < i log |l + ^ ^'*»^^| ^ ^ 



ies ieS 



A typical capacity region of Gaussian MAC for 3 users is illustrated in Figured 

As in [14 , we will unify our treatment and present the results for MAC in more general terms. For 
any subset S* C and r = (ri, . . . , r^), we use the shorthand notation 



r{S) := E 

for the summation of components of r with indices in S 



ies 
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For any function g mapping subsets of Q to R+, let V{g) denote the polyhedron 



{r e : r{S) < g{S), VS* C 



If the function g satisfies (i) g(0) = 0, (ii) g{S) < g{T) if 5 C T, and (iii) g{S)+g{T) > g{Sr]T)+g{SUT) 
for all subsets 5* and T of fl, then g is called a rank function. Property (ii) and (iii) are called the monotonia 
and suhmodular property respectively. We say that the polyhedron V{g) is a polymatroid when 5 is a rank 
function. 

The capacity regions of scalar and vector Gaussian MAC are polymatroids. For scalar Gaussian MAC, 
we define the rank function as 

5(^) = ilog(l + l^P.). (1) 

ies 

On the vector case, we define the rank function as 

ff(^) = ilog I+l^P,s,sf . (2) 



N 



The collection of points that achieve equality on total data rate r(n) = g{^l) is called the dominant 
face. Given a vector r in V{g), we say that the set S* is a bottelneck of r if 

riS)=giS). 

For any function g : 2^^ — > M_|_, if 

g{A) < g{B) implies g{A U C) < g{B U C) 

whenever Ar\C = ^ = B nC and \A\ = \B\, then we say that g satisfies the order property. The heuristic 
meaning is as follow. If the achievable sum rate of group A is less than or equal to that of group B, then 
even if they cooperate with users in group C, the sum rate of group A U C is still less than that of the 
group B UC. 



Example 1 Consider the function g defined as g{{a}) — 1, g{{b}) — 2, g{{c}) ~ g{{b, c}) — g{{a, b}) = 3, 
g{{a,c}) = g{{a,b,c}) = 4. We have g{{a}) < <?({&}) but g{{a,c}) > g{{b,c}). This function does not 
satisfy the order property. 

The scalar Gaussian MAC satisfies the order property, but the vector Gaussian MAC in general does 
not. 

In both the scalar and vector Gaussian MAC, there is a canonical choice of disagreement point for the 
Nash bargaining solution. Each user can treat the signal of the others as noise and decode independently. 
For scalar Gaussian MAC, the resulting data rate for user i is 

In the vector case, if user i uses linear MMSE receiver with no joint processing with others, the data rate 
is [15] 

c(^P,sfM-is,), 

where is the matrix 

1+ — y^P,sfs,. 
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In both cases, they can be expressed in terms of the rank function as 

d*:^gin)^gin\{t}). (3) 

The rate vector d* ~ . . . , d*j^) is called the canonical disagreement point. 

The next lemma is a useful consequence of the order property. The proof is straightforward and is 
omitted. 

Lemma 5. Let g : 2^ — > be a function that satisfies the order property. 

1. If g{{\}) < g{{2}) < ... < g{{K}) are in nondecreasing order, then g{{l,2, . . . ,i}) < g{S) for all 
S' C o/ size i . 

2. For any subset A, the function g'{S) — g{S U A) — g{A) defined for S C n\A also satisfies the order 
property. 

4.1 Symmetric Rate Allocation 

The computation of the symmetric capacity in V{g) amounts to finding the tightest constraint among 
'''{S) < g{S) for all subsets S. Each component of the symmetric rate cannot exceed g{S)/\S\ for all 
S C il. The symmetric rate allocation can be computed by checking 2^ — 1 constraints. 

If g satisfies the order property, the symmetric capacity can be computed more efficiently using 
Lemma [H The computation requires taking the minimum of only K numbers. 

Theorem 6. Let g : 2^^ — > K+ be a function that satisfies the order property. By relabeling we can assume 
that g{{i}) < gi{j}) whenever i < j. The symmetric capacity in the polyhedron V{g) equals 

K-mm[p{{l,...,k}): = 

Corollary 7. The infimum of fair efficiency of symmetric fairness, taken over all K-user MAC, is zero. 
Proof. Suppose that the power of user 1, Pi, is much less than the others, so that 

.g({l}) = min{i5({l,...,fc}): fc = l,...,i^}. 

We see that the symmetric capacity approaches zero when Pi approaches zero, while the maximal sum 
approaches a positive constant. □ 

We next compare two MACs with different power constraints. If the power constraints Pi, . . . ,Pk 
become more disperse, then the symmetric capacity will decrease. 

Theorem 8. Let Csym^P ^ N) be the symmetric capacity of a scalar Gaussian MAC with power constraints 
P and noise power N. If F ^ P', then Csym(P, N) > Csym{P',N), i.e. the symmetric capacity of 
multiple- access channel with fixed noise power is a Schur-concave function. 

Proof. Without loss of generality, we assume that the power constraints are sorted in nondecreasing order. 
Pi < P2 < ■ ■ ■ < Pk , and P^ < P2 < . • . < P^ . Since scalar Gaussian MAC satisfies the order property, 
the corresponding symmetric rate allocations r*^™ and s*^^™ are given by 

rf™=min{ic(l^PO:fc = l,...,if} 

and 

.f™=min{ic7(l^/^')^fc-l>---^^} 
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for all i. As 

k k 



i=l i=l 



for all k, we have 



A; 
i=l 

Therefore rf^™ > sf^'". □ 



4.2 Max-min and Proportional Fair Rate Allocation 

The following is a useful characterization of max-min fairness |16) . 

Lemma 9. For any function vector r is max-min fair in the polyhedron V{g) if and only 

if for all i, the ith component is largest in some bottleneck. In other words, for i = 1, . . . , K , i is contained 
in a bottleneck B and ri > rj for all j E B . 

The proof of Lemma [9] is contained in the Appendix. The next theorem is the main theorem in this 
section. The proof of Theorem [TO] is in the Appendix (See Prop. [55] and [551) 

Theorem 10. Suppose that g : 2^^ — > M_|_ satisfies the submodular property. The max-min fair point r*^*^ 
in 'P{g) is on the dominant face, and is majorized by every point on the dominant face. 

It is noted that in the theorem we do not assume that the function g is a rank function. The result 
holds as long as we have the submodular property. The max-min fair allocation is fairer than any other 
point on the dominant face in the sense of fairness induced by majorization. The max-min fair solution 
also has the following interpretations. 

Corollary 11. If g : 2^ ^ TZ^^ satisfies the submodular property, the order property and gifli) — 0, then 
the max-min fair (and hence the proportional fair) solution r* maximizes r[x], ^[x] +^[2]: ■ ■ ■ ^ o.'^d '"'[j] 
simultaneously over all points on the dominant face. 

Corollary 12. If g satisfies the submodular property, then the max-min point and the proportional fair 
point in V{g) coincide. 

Proof. Since log is a concave function, ^j^^og{ri) is Schur-concave by Lemma [TJ Because the max-min 
point r*^*^ is majorized by any point r on the dominant face, we have 

^log(rf"^)>5:iog(rO. 

i i 

Hence r^^^r™. □ 

Corollary 13. If g satisfies the submodular property, then the max-min (and the proportional fair) point 
in 'P{g) is the point on the dominant face that minimizes the Euclidean norm. 

Proof. The function f{x) — x^ is convex. The proof is similar to the proof of the last corollary. □ 

Corollary 14. // g satisfies the submodular property, the fairness efficiency of max-min fairness and 
proportional fairness is equal to 1. 

Theorem 15. Suppose that g : 2^^ ^ K+ satisfies both submodular and order property, and assume 
^({l}) < 5({2}) < ■■• < 5({^}) after suitable relabeling. Then the components of max-min fair (and 
hence the proportional fair) solution r* in T'{g) are in nondecreasing order, r\ < < • . • < • 
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Proof. Consider any i in fl. The index i is contained in a bottleneck Ai of r* so that r* — max{r* : j e Ai}. 
If i — 1 G Ai, then r*_^ < r* . Otherwise, suppose i — 1 ^ Ai and let S denote Ai \ {i}. Then 

r*{S U{i- 1}) < g{S U {i - 1}) < g{S U {i}) = r*{S U {i}). 

We have used the order property in the second inequality. This implies that 7'*_j^ ^ CH 

A typical class of functions that satisfies the order property is the generalized symmetric functions. A 
rank function g is said to be generalized symmetric if it has the form 

giS) - miS)), (4) 

where </> is a monotonic increasing and concave function with (f>{0) — and Q G R:^. The rank function 
in the scalar MAC is an example of generalized symmetric function. 

The next theorem compares two Gaussian MACs, with the same total power but different distribution 
in the power constraints. It shows that if the distribution of power constraints is more spread out, so does 
the corresponding max-min fair rate allocation. The proof is relegated to the appendix. 

Theorem 16. Let g be a generalized symmetric rank function defined as g{S) = (f>(P{S)), for some 
vector P G K;^. Let P be a vector that majorizes P, and g' be the generalized symmetric rank function 
g'{S) — (t>(P{S)). Then the max-min fair capacity associated to P is the same as the max-min fair capacity 
associated with P, and he max-min fair point in 'P{g) is majorized by the max-min fair point in V{g'). 

Corollary 17. Let r*(P,7i) be max-min fair in a Gaussian MAC with power constraints Pi < . ■ . < Pk 
and noise power n. LfP ^ P', then r*(P,n) ^ r*(P',n). 



4.3 Algorithm 

We present a general recursive algorithm that computes the fair solution in Gaussian MAC. It is a variation 
of the algorithm in ^16, p. 527], which computes the max-min fair rate vector in flow control problem. The 
algorithm to be described below exploits the submodular property, and has shorter running time. The 
basic idea is contained in the next proposition. 

Proposition 18. Let g be a function mapping 2^ to R-|_, and let Sq be a subset of ft that achieves the 
minimum 

min q(S)/\S\. 

Let r* denote the the max-min or proportional fair point (they are the same by Corollaru \lS\) in 'P{g), we 
have r* > g{So)/\SQ\ for all i, with equality when i E Sq. 

Proof. For each i, i is contained in a bottleneck Bi of r*, so that > rj for all j G -B,;. So must be 
larger than or equal to the average r*(_B,j)/|_B,j|, and thereby 



B. 

Therefore r* > g(S'o)/|>S'o| for all i. Summing over all i G 5*0, we obtain 

r*{So) > g{So). 

We must have equality in all the above inequalities. In particular, r* — g(>5'o)/|S'o| for all i G Sq. The set 
So is in fact a bottleneck of r* . □ 

This motivates the max-min algorithm. 
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M£LK-min algorithm The algorithm starts by first obtaining the subset So £ Q described in Proposi- 
tion[THl and set the rate of users in 5*0 to g{So)/\So\. The rate of other users are computed by recursively 
applying the above computation to fl' := ^ \ Sq with 

g' {S):^g{S Li So) -g{So) (5) 

for S c n'. 

Remark: It is noted that we only need the submodular property in proving the correctness of the 
max-min algorithm. The function g need not satisfy the monotonic property or g{0) = 0. We will use the 
following lemma in proving the correctness of the algorithm. Note that the lemma holds in general for 
arbitrary g. 

Lemma 19. Let g be a function from 2^ to M_|_, and So be chosen such that 

giSo)/\So\ = ^min^g(5)/|S'|. 

Define the function 

g\S) -.^ giS U So) - g{So) 

for sen' ■=n\So. 

1. (Non-negativity) g' > 0, for all S C fl' . 

2. (Extension of bottleneck) Let r be a vector in T'{g) such that r(S'o) = g{So)- Let r' be the restriction 
of the vector r on fl' . If B' is a bottleneck of r' in Vig'), then B' U So is a bottleneck of r in 'P{g). 

3. (Preservation of order property) If g satisfies the order property, so does g' . 
Proof. (HI) By construction, we have 

\SUSo\ 



l(.SuSo)>giSo)^ 



\So\ 



Hence g{S U So) > g{So). 

(m r'{B') equals g'{B'), which is g{B' U So) - giSo) by definition. Therefore 

g{B'liSo) = r'{B')+g{So) 
= r(S')+r(5o) 
= r{B' U So), 

and B' U j4 is a bottleneck of r in P{g)- 

© Let S, T and U be subsets oin\So so that l^] = |T| and 5 n ?7 = = T n f/. If g'{S) < g'(T), 
then 

g'{SUU)^g{SUUUA)-g(A) 
< g{T UUUA)- g{A) 
^g'{TUU). 

□ 

Theorem 20. Suppose g : 2^ — » IR+ satisfies the submodular property. In the polyhedron V{g), the result 
obtained by the max-min algorithm is the max-min fair vector. 
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Proof. Let r be the vector returned by the max-min algorithm, and let 5*0 be a subset of 17 such that 
g{So)/\So\ = min0-^5csi 5'('S')/|'S'|. If 5*0 = fl, then the components in r are constant and equal g{So)/\So\. 
It belongs to P{g) because for any S C il 

riS)^\S\^<g{S). (6) 

It is easy to see that r is max-min fair. 

Otherwise, if 5*0 C $1, then we have to apply the algorithm recursively. In this case, the vector r 
satisfies the following properties: (i) = g{SQ)/\So\ for all i (z Sq and (ii) the components of r with index 
in f2 \ 5*0 yield the max-min solution to the polymatroid M.' on fl\ So with rank function g' defined as 
in (HI. It is easy to check that g' also satisfies the submodular property. We first verify that r is in 7^(g). 
For any subset A C J7, we can write A a.s AiU A2 with Ai C 5*0 and A2 fl 5*0 = 0. We decompose r{A) as 

r(A) =r(Ai)+r(A2) 

<\A^\^+giA2USo)~g{So) 
<g{Ai)+g{A2USo)-giSo) 

In the last inequality, we have used the defining property of So, i.e., |v4i|g(S'o)/|S'o| < g{Ai). By the 
submodularity of we have 

5(5-0) + giAi U A2) > g{Ai) + g{A2 U ^o). 

Therefore r{A) < g{Ai U A2), and thus r is in V{g). 

We now show that for each i = 1, . . . , if, i is in some bottleneck Ai such that = max{rj : j e Ai}. 
We will apply Lemma[9]and conclude that r is the max-min vector in P{g). For i G So, we can take So 
as the required bottleneck Ai. For i ^ So, i is an element of some bottleneck B' in the polyhedron P(.g') 
such that ri > rj for all j € B' . By part ^ in the previous lemma, B' U 6*0 is a bottleneck of V{g'). By 
Prop.[THl we can show that > rj for all j & So. Indeed, 

^ 9'iB') 

1 / , \SoUB'\ . 

_ gjSo) _ 

Therefore i is in the bottleneck 5*0 U B' and > rj for all j £ B'USq. The vector r is thereby max-min 
fair by Lemma [9] □ 

Example 1 (continued) We compute the max-min fair vector in Example 1. The minimum 

min q(S)/\S\ 

is achieved when S = {a}. We set rf^*^ — g{{a]) — 1. Next define 

g\S) ■.= g{{a}yjS)~g{{a}) 

for S C {6, c}. Now g'{{h}) = 2 and g'{{c}) = 3 = g'{{h,c\). So in this recursive step, we have the 
minimum 

5'({6,c})/2-3/2. 
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The resulting max-min fair solution is 

^MM ^ (1^3/2^3/2). 



In the max-min algorithm, if we compute the minimum of g[S)/\S\^ 7^ C 17, in a straightforward 
manner by comparing 2^ — 1 numbers, the complexity of the algorithm is exponential in the number of 
users. A more efficient implementation was described in |^ if g is a rank function. However, when the 
function g satisfies the order property, we have a much faster algorithm. 

Proposition 21. // the function g satisfies the order property and the submodular property, then the 
max-min fair point in V{g) can be computed in 0{K^) time. 

Proof. In the max-min algorithm, instead of finding the minimum 

min q{S)l\S\ 

over all subsets of fi, we sort g{{l}),g{{2}), . . . ,g{{K}) in nondecreasing order. This can be done in 
0{K\og{K)) time. For notational convenience, we relabel the users so that 

g{{l])<g{{2])<...<g{{K}). 

Since function g satisfies the order property, for any k, the minimum 

min{5(S')/fc: 5 C 17, IS*] = fc} 

is achieved by {l,2,...,fc} by Lemma [5] Instead of comparing g{S)/\S\ over all subsets S of il, it is 
sufficient to examine (7({1, . . . , fc})/fc, for k — 1, . . . ^K. The minimum can be found in 0{K) time. In the 
next recursion, the function g' {S) — g{S U Sq) — 5(5*0) also satisfies the order property, and 

.9'({fco + 1}) < g'i{ko + 2}) < . . . < g'{{K}). 

The recursion can continue without any further sorting. There are at most K recursive steps and each 
step takes 0{K) time. The total complexity is therefore 0{K^). □ 



Example 2 Consider a scalar Gaussian MAC with 4 users. Their power constraints are 2, 8, 200, and 
300. We let P be the vector (2,8,200,300). The noise power at the receiver is equal to 1. Let 

5i(5) :=0.51og(l + P(5)) 

for 5* C J7. Here we use the natural logarithm function. We want to find the max-min fair rate allocation 
or the proportional fair rate allocation in 'P{gi). The function gi is a rank function and satisfies the order 
property. 

We first compute the minimum of .gi({l}, gi({l,2})/2, 5i({l, 2, 3})/3 and gi({l, 2, 3, 4})/4. The min- 
imum is .91 ({1}) = 0.5493. We set = 0.5493. 
In the next recursive step, set 

32(5) :=5i(5U{l})-5i({l}) 

for S C {2,3,4}. The minimum of g2({2}), g2({2,3})/2 and g3({2, 3, 4})/3 is g2({2}). We set rf*^ = 
52({2}) = 0.6496. 
Let 

<?3(5) 52(5 U {2}) - g2({2}) = gi{S U {1, 2}) - 5i({l, 2}) 
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for S C {3,4}, and compute the minimum of g3({3}) and (;3({3, 4})/2. The minimum is (73({3,4})2 = 
0.9596, and we assign 0.9596 to both rf"^ and r|"^. 
The max-min fair rate aUocation is thus 

J.MM ^ (0.5493,0.6496,0.9596,0.9596). 

The computation of the Nash bargaining solution in 'P(g) with disagreement point d amounts to finding 
the max-min fair solution in V^g'), where 

g'{S)^giS)-d{S). 

It is noted that if 5 is a rank function, the translated g' in general does not satisfy the monotonic property. 
However, the max-min algorithm works without assuming the monotonic property. We can apply the 
max-min algorithm to find the Nash bargaining solution for any disagreement point. 

We conclude this section by presenting an algorithm for computing Nash bargaining solution when the 
rank function is generalized symmetric. 

Lemma 22. Let g be a generalized symmetric rank function on fl, 

g(5)-0(P(5)), 
for sen, and P is a vector in R;^ such that. 

Pl<P2< ■■■<Pk. 

Then for any fc G {1, . . . , K}, the minimum 

min |.g(S') 4- ^ 5(1^ \ {i}) : 5 C 17, |5| = fc} 

is equal to 

k 

g{{l,2,...,k}) + J2gin\{z}). 
1=1 

Proof. Let Pi denote X^j^i where the summation is over all indices except i. The lemma claims that 
(f>(P{{l, . . . , k}) + X]i=i 'P{Pi)j is the minimum. By Lemma[T] it suffices to show that 

{PuP2,...,Pk,Pi + ... + Pk) 

is majorized by 

{Pii 7 Pi2 7 • • • 5 Pik ' Pil + ■ • • + Pik ) 

for any choice of ii, . . . ,ik, 1 < ii < ■ ■ . < ik < K . 

By subtracting J^iLi Pi from both vectors, we only need to show that 

{PuP2,---,Pk,Pk+l + -.- + PK)h [P^^ , P., , . . . , P,;, , ^ P, 

where B denote the set {ii, . . . , ik}- 

Let Q be the vector ^P^^ , P^^ , . . . , Pi^^ , J2ji^B ^i) ' ^^"^ Q[2], ■ ■ ■ , Q[k] be the components of Q in 
nondecreasing order. It is easy to see that 

k 

i=l 

Therefore (Pi, P2, . . . , Pk,Pk+i + . . . + Pk) majorizes Q. This finishes the proof of the lemma. □ 
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Proposition 23. For a generalized symmetric rank function g, the Nash bargaining solution in V{g) with 
the canonical disagreement point can he computed in 0{K'^) time. 

Proof. Assume without loss of generality that the power constraints are arranged in nondecreasing order. 
Let d* denote the canonical disagreement point, i.e., iov i = \, . . . ,K, dl := g{il) — g{il, \ {i}). Let 

hiS) := g{S) - d*{S) = g{S) + ^ g(n \ {i}) - \S\g{n). 

ies 

We relabel the users so that h{{l}) < . . . < h{{K}). By the previous lemma, the minimum 

min{/i(S'): 5 C l^j = fc} 

is h{{l, . . . , k}). Therefore, 

mm h(S)/\S\= mm h({l, . . . ,k})/k. 
0#sca ^''11 k=i,...,K 

The minimum can be obtained efficiently after sorting /i({l}), . . . , h{{K^). Suppose that the minimum is 

/i({l,...,io})/|io|. Wesetrf^ to • • • ,'»o})/No| + < for i = 1, . . . , io- 

We next show that the same procedure can be repeated in the next recursive step. Let Sq = {1, . . . , io}, 
and g'{S) = g{SLiSo)—g{So) for 5* C fi' := {io+1, • • • , K}. It is easy to see that <?' is generalized symmetric. 
Also, we can verify that • • • , ) canonical disagreement point for g'. Indeed, 

d*=gin)-g{n\{z}) 
= g'm-g\n'\{z}) 

for all i ^ Sq. 

Each recursive step takes 0{K)). As there are at most K steps, the complexity for computing the 
Nash bargaining solution with the canonical disagreement point is 0{K^). □ 

Example 2 (cont'd) We compute the Nash bargaining solution in the MAC as in Example 2, with the 
canonical disagreement point d*. 

dl = gi({l, 2, 3, 4}) - gi({2, 3, 4}) = 0.0020 
d*2 = 5i({l, 2, 3, 4}) - gi({l, 3, 4}) = 0.0079 
dl = 5i({l, 2, 3, 4}) - gi({l, 2, 4}) = 0.2483 
dl = 3i({l, 2, 3, 4}) - <?i({l, 2, 3}) = 0.4423 

The minimum of 

{g^{{l,...,k})-dl-...-dl)/k 
for k = 1,2,3,4, is ffi({l}) — dl. Therefore, 

rf"^ = {gii{l})- dl) + dl = 0.5^93. 

For S C {2,3,4}, let g2{S) := gi{S U {1}) - ffi({l}). 

52({2}) ~d*2= 0.6418 
g2{{2,3}) - d*2 - dl 



2 

g2{{2,3A})-d*2-dl-dl 



0.9352 
0.6235 
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The last equation yields the bottleneck. We set 

^MM ^ Q g235 + rf* = 0.6314 
^MM ^ Q g235 + rf* = 0.8718 
^MM ^ Q g235 + rf* = 1.0657 

The resulting Nash bargaining solution is 

(0.5493, 0.6314, 0.8718, 1.0657). 

5 Broadcast Channels 

In a -ff-user Gaussian broadcast channel, the received signal of the ith user is 

Yi=X + Zi 

where X is a zero-mean Gaussian random variable with variance Pt and Zi is the noise at the ith 
receiver, which is modeled as a Gaussian variable with mean zero and variance Ni. We will assume that 
Ni < N2 < ■ . ■ < Nk throughout this section. Every point r on the boundary of the capacity region 
satisfies 

'^j = o log H ^^r^i ) 

for some ai , . . . , ax such that 

K 

^ttj = 1, Uj > 0. 

5.1 Symmetric Capacity 

In order to obtain the symmetric capacity in a BC, we solve a related problem of finding the power 
distribution so that the users have a common SINR 7, i.e., 

_ P]_ _ P2 _ _ PK_ 

^~ Ni ~ N2+P1 Nk+Pk-1 + ...+P1' 

The noise vector N can be expressed in terms of the power vector p by a matrix multiplication 







-1/7 


" 






N2 




-1 1/7 







P2 


Nk. 




.-1 -1 


-1 1/7. 







Let G-y denote the lower triangular matrix in the above equation. The diagonal elements of G-y all equal 
1/7, and the elements below the diagonal are all —1. The following lemma is obtained by straightforward 
calculation. 

Lemma 24. The inverse of is a non-negative matrix. The {i,j)-entry of G~^ is 

ifi<j, 

7 « = j, 

7^(1 + 7)*^-'^"^ otherwise. 
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Hence, given the noise powers N = (A^i, . . . , Nk) and the SINR requirement 7, we get the corresponding 
power aUocation by multiplying by N. The following theorem is an immediate consequence. It says 
that in a Gaussian BC, users with lower noise power uses less power. 

Theorem 25. In a Gaussian BC with noise powers Ni < . . . < Nk, the powers corresponding to the 
symmetric rate allocation are in increasing order, 

sym ^ sym ^ ^ sym 

Pi <P2 < ■■■<Pk ■ 



Proof. Let p be the power vector so that all users have SINR 7, and let h^j be the (z,j)-entry of G^^. 
The power of user i can be obtained by 

Pi = hiiNi + . . . + huNi. 
Since the first term is nonnegative, we can remove the first term and get 

> h^2N2 + haNs + ... + huN. 
It is clear from the previous lemma that hij depends only on i — 7, hence 

Pi > /''(i-l)lA^2 + ^(i-l)2^3 + • . • + /l(i-l)(i-l)iVi 

We now use the assumption A^i < . . . < Nk to obtain 

Pi > ^(j^l)lA^l + ^(j-l)2^^2 + • • . + /l(,_i)(j_i)iVj-l 

= Pi-i- 

Consequently pi < P2 < • ■ • < Pk- CH 
We denote the required total power by 0(N, 7), which can be computed by 

^(N, 7) = [!,..., 1].G;1N 

= NiOii-,) + A^202(7) + • • • + NkOk^) 

where ^^(7) denote the sum of elements in the jth column of G^^. It is noted that the function ^^ (7) is a 
convex function of 7 for all j. Furthermore, we have 

^i(7)>^2(7)>--->^k(7)- (7) 

Given a noise vector N, the function is a convex and monotonically increasing function of 7. 

The next theorem compares the symmetric capacity of two broadcast channels with the same total 
power constraint. 

Theorem 26. // N ^ N', then Csj^™(N, Pt) < CsymC^' , Pt) ■ In other words, Gsy,n{-,PT) is a Schur- 
convex function. 

Proof. Assume without loss of generality that the components of Ni and N2 are sorted in nondecreasing 
order. For a given total power constraint and noise vector N, we can obtain the SINR by solving the 
equation 

0(N,7) =Pr. 

The proof is complete if we can show that (/)(N,7) > (/)(N',7) for all 7, i.e., </)(•, 7) is Schur-concave for 
all 7 (Fig. [3]). Indeed, for i > j, we have Ni > Nj and 

by ([7]). Therefore (f>{-,"f) is Schur-concave by Schur's criterion. □ 
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Figure 3: The total power required for symmetric rate allocation as a function of SINR. The noise vector 
associated to the lower curve majorizes the noise vector associated to the upper one. The dash lines show 
how to find the SINR when the total power is given. 

Theorem 27. For fixed noise powers, the symmetric capacity is a convex function of the total power. 

Proof. For a fixed noise vector N, the function (/)(N,7) is a convex function of 7. Hence the inverse 
function is a concave function. □ 

Algorithm We have a numerical algorithm computing the symmetric capacity in Gaussian BC by means 
of the function (/)(N,7). For a given total power Pt, we search for the value of 7 so that (/)(N,7) = Pt. 
This can be done easily as is a monotonic function of 7. We then compute the data rate from 7. 

Corollary 28. The infimum of r] gym, taken over all K-user Gaussian BC, is zero. 

Proof. The sum capacity is attained if we allocate all power to user 1, 



Suppose the noise power of user K, Nk, is increased, while the others arc fixed, the value of the fimction 
(/)(N, 7) is increased for all 7. Then the symmetric capacity is decreased, but the sum capacity remains 



5.2 Proportional Fair Capacity 

The capacity region written in the following form is useful for computing the proportional fair allocation: 




constant. By taking Nk approaching infinity, rjsym approaches zero. 



□ 




for i = 1, . . . , ii', 



(8) 



with 



Q = Xo < Xi < X2 < ■ ■ ■ < XK-1 < Xk = 1- 
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The quantity XiPr represents the sum of powers Pi + . . . + Pi. 

Theorem 29. Let "I>i he a concave and monotonically increasing Junction for i — 1, . . . ,K. Using the 
notation in fS]), the point that maximizes 

K 

/(r) :=5]<i>.(r,) 

i=l 

in the capacity region of a BC with total power Pt and noise powers Ni < . . . < Nk satisfies 

(9) 



Ni + x,Pt N,+i + x^Pt 

for I ^ I,..., K -I. 

In particular a proportional fair allocation satisfies the equation 

r,{N, + x,Pt) = r,+i(iV,+i + x^Pt) (10) 

fori = l,...,K-l. 
Proof. Suppose that 

V2 ^^N, + x,^iPt' J i=i.....,K ^ ' 

maximizes /(r). This point hcs on the boundary of the capacity region. The tangent plane at r must be 
orthogonal to the gradient 

V/(r)= (<i>;(ri),$^(r2),...,$i,(rfe)) 

i.e., it must be orthogonal to r' — r for all r' on the tangent plane. 
Differentiate pip with respect to a;^, we obtain for i = l,...,i4r — 1, 



Prp Prp 



Ni+XiPr' N,+i+x,Pt 



The two fractions in the above vector is in the ith and (i + l)st component. This must be orthogonal to 
the gradient V/(r). Equating the dot product to zero, we get 

1 I P-r 

2 Ni + XiPr 2 Ni+i + XiPr 

For the proportional fair point, we take <I>i to be the log function for all i. □ 

Corollary 30. Let r^^ and be the rate vector and SINR vector corresponding to the proportional 

fair allocation, with Ni < . . . < N^. Then the rates and SINR of the users are in decreasing order, i.e., 

rPF > ^PF > > PF 

and 

7r>72^^>...>7r. 



Proof. From pTH) . we have 



N,+i + x ,Pt ^ 
Ni + XiPr 



Since the rate is a monotonically increasing function of the SINR, the inequalities about SINR follows 
immediately. □ 
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Corollary 31. The infimum of rjpF, taken over all K-user Gaussian BC, is equal to the lower bound 
1/K. 

Proof. Suppose that we fix the noise powers Ni, . . . , Nx and take Pt — > 0. Equation (fTU|) imphes that 

rrN. « rf^N, 

for all i and j. When Pt is small, 

rP^ = i log f 1 + ^ 1 « 

Hence pf^ « p^"^ for all i and j, and pf^ Pt/K at the proportional fair allocation as Pt 0. We 
have the following limits, 



r - — 



and 

CpF 



K 



1 " PF 

l — l 

~ V — 

We obtain 

Cpf 1 / iVi N,^ 

The right hand side can be arbitrarily close to 1/K ii Ni <C A'i for all i = 2, 3 . . . , if. □ 

Theorem 32. In a Gaussian BG with noise power Ni < . . . < Nk, the corresponding powers correspond- 
ing to the proportional fair allocation are in increasing order, 

PF PF PF 

Pi <P2 < ■ ■ ■< Pk ■ 



Proof. We first prove the theorem in a two-user case. The noise power of user i is Ni, {i = 1,2) with 
A^i < N2. By Theorem \29\ the proportional fair rate allocation is 

1, fNi + aPT\ 
"^=2^°H^V^j 

1 , ( N2+Pt \ 

''^^2^°H^^^— ^j, 

where a, < a < 1 is chosen such that 

ri{Ni + uPt) = r2{N2 + uPt). 

That is, the value of a satisfies 

[Ni + aPT) log [^^^^^) - {N2 + aPT) log (^|^) (12) 

The power of user 1 is uPt and the power of user 2 is (1 — a)PT. We want to show that (fT2l) cannot hold 
if a > 0.5. 
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We will apply the inequality 

a; log (l + -) < ^ 

which holds for x > maxjO, —6}. This inequality is an immediate consequence of the inequality 

< 1 + x, Va; 7^ 0. 
Applying the inequality, we get an upper bound 

R.H.S. of (ini) = (iV2 + aPr) log f 1 + ^}~'^^^^ ) < (1 - a)PT, 



N2 + aPr 



and a lower bound 



L.H.S. of (HH) = -{Ni + aPr) log (1 + — i^^^ ) > -{-aPr) = aPr- 

\ Ni+ aPr ) 

If a > 0.5, we can combine the two bounds above, 

/ iVi + aPr \ 
{Ni + aPr) log ( — -j > aPr > (1 - a)PT 

Therefore, we must have strict inequality in (|12p when a > 0.5. It is noted that the bounds in the above 
argument is valid for any A^i and N2- 

In a Gaussian BC with K users, consider the pair of consecutive users i and z + 1. At the proportional 
fair point in the capacity region. 



P^+l 



by Theorem [29l By setting 



we get 



Ni := TV, + E 

i=i 

j-i 

Pt := P, + P.+i 
a:=P,/(P, + P,+i), 



[ Ni+ aPr \ , ~ ~ , [ N2 + Pt 
{Ni + &Pt) log — ^-^ \^{N2+ &Pt) log ' 



Ni I \N2 + aPr I 

We can proceed as in the 2-user case and conclude that pi < pi+i- □ 

From numerical examples, we have the following conjecture. 

Conjecture 33. Let Cpi?(N, Pt) denote the proportional fair capacity of a Gaussian BC with noise vector 
N and total power Pt- If N ^ N', then Cpf{^,Pt) < CpfCN^Pt) 
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5.3 Algorithm 



We will use the notation for the capacity region in ([5]) , and present an numerical algorithm that maximizes 
Si=i ^i(^i) in the capacity region of a Gaussian BC. It is assumed that is a strictly monotonically 
increasing and concave function, so that the inverse of the derivative <i>^ is easy to compute. If is the 
logarithm function, then the resulting point is the proportional fair solution. If $i(a;) = log(a; — di), the 
result is the Nash bargaining solution with disagreement point d. 
Theorem says that the equation 



$^+i(r,+i) = $^(r,) 



N., + x,Pt 



must holds for the optimal solution. We use this equation to express x^+i in terms of Xi for i — 2, . . . , K — 1. 
It reduces the problem to a one dimensional search. Given any xi we first compute r2 by 



N2 + XIPT 

' Ni+xiPt 
and get X2 by solving 

1, (N2 + X2PT 
^2 = 7: log 



2 °\N2+xiPt- 

In similar way, we compute and Xi for i = 3, 4, . . . , if — 1. Finally rx can be obtained once we know 

XK-l- 

Define the function xi^i) 

X{xi) := <^'ii{rK){NK~i + xk-iPt) - ^'k-i{i-k-i){Nk + xk-iPt) 

It is a function of xi as the variable tk, tk-i, xk and xk-i all depend on xi. We can now search for 
the zero of xi^i) numerically, say xi^i) = 0. From x\^ we get x* by method described above. The vector 
(x*, . . . tX*j^) will satisfy the condition in Theorem 1291 and hence is the optimal solution. Since any zeros 
of X gives rise to an optimal solution and we know that there the optimal solution is unique, the function 
X has only one zero. 



6 Conclusion 

We show how to pick a point in the capacity region of Gaussian MAC and BC according to some fairness 
criteria. In the Gaussian MAC, there is a strong notion of fairness, namely there is a point on the dominant 
face that are majorized by all other points on the dominant face, and are both max-min and proportional 
fair. We can thus call this the fair point in the capacity region. In some particular cases, the fair point can 
be computed in 0{K^) time. For the Gaussian BC, the problem of locating the proportional fair solution 
or Nash bargaining solution reduces to a one-dimensional search. In both channels, fair rate allocation 
can be compute efficiently. 



7 Appendix 

Proof of Lemma\^ (<^=) For any f G il, suppose that i is contained in a bottleneck B, and — max{rj : 
j S i?}. If we want to increase r^, we have to decrease for some other k Cz B. Since r.i is the largest in 
{rj : j £ i?}, we must have < r^. The vector r is thus max-min fair. 

(=>) Conversely, suppose that r is a vector such that Bi, . . . , Bt are all the bottlenecks that contain 
«, and Ti is not maximal in all such bottlenecks, i.e., < max{rj : j £ Bt] for alH = 1, . . . ,T. We can 
choose it in Bt such that r^^ > r^. If we increase by e and decrease each ri^ by sufficiently small e, the 
resulting vector remains in V{g). The vector r is thereby not max-min fair. □ 
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Lemma 34. Suppose that a function g : 2^ ^ satisfies the submodular property. Union and intersec- 
tion of two bottlenecks of r in 'P{g) are also bottlenecks of v. 

Proof. Suppose that S and T are both bottlenecks of r, i.e., r{S) — g{S) and r(T) = 5(5'). 

r(5)+r(r)-g(5)+5(T) 

>g{SUT)+giSnT) 
>r(5UT) + r(5nT) 
= r(5)+r(r) 

Therefore, aU inequahties above are in fact equahties. In particular, r(5' LiT) — g{S U T) and r{S PiT) — 
giSnT). □ 

The proof of Theorem 1101 is divided into the next two propositions. 

Proposition 35. Suppose thatr'^^^^ be the max-min fair point inV{g), where g satisfies the order property. 
By relabeling, we can assume without loss of generality that 

MM _ _ „MM MM _ _„MM ^ ^MM _ _ „MM o^ 

where 

< ii < i2 < ■ . ■ < iL = K. 
The sets {!,..., i^}, £ — 1, . . . ,L, are bottlenecks of r*^*^ . In particular, we have 

i.e., the max-min fair solution lies on the dominant face of V{g). 

Proof. Let Bi denote the set {1, . . . , } for £ = 1, . . . , L. 

For any element j G Bi, there is a bottleneck Aj so that rj = maxfr^ : i G Aj}. If j S Bi and k ^ Bi, 
then rj'^^'^ < r^^ . Hence Aj must be a subset of Bi, for all j £ Bi. By taking the union of Aj over all 
j eBi, we get 

^1 = U ' 

and we can conclude that Bi is also a bottleneck of r*^*^ by Lemma [Ml 

By similar argument, we can show that B( is bottleneck of r^^^^ for all £ £ {!,..., L}. □ 

Proposition 36. For £ — 1, . . . , L, let B^ be the set {1,2,..., ig}, with < ii < i2 < . . . < ii = K. Let 
w be a vector in M.^ such that 

Wi = . . . = Wii < Wii+i = . . . = < . . . < = ...^Wi^ 

All vectors in the region 

7e {r £ : r(B,) < w(B,), £^1,...,L} 

satisfy 

k k 

j=i i=i 

for k — 1, . . . , K . Consequently, w is majorized by all points in TZ such that r(f2) = w(f2) 
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Proof. Suppose that v is a point in TZ that does not majorize w. There is an index /c, 1 < fc < if , such 
that 

k 

Wi 

We can find an i so that ig < k < ig+i. (Define ia := and Bq := if necessary.) 

Consider the collection T of all subsets S C ft such that Bg C S C Bg^i and \S\ = k. The number of 
such subsets is 

\ k - ig 

Since Wi is constant for i G i?£+i \ -Bf, we have 
for all S eT. Hence 



k 

' Wi 



We sum the above over all S 



Wi 

ies i=i i=i ieS 



EE«^>EE"^- (14) 

The left hand side in the above inequality equals 

a^^v^ + b ^ Vi, 

i—l i—i^+l 

where 

^ aik - ii) 
ii+i - ii 

Similarly, the right hand side of ([T4| equals 

if if+l 

+ 6 ^ Wi. 

We rewrite as 

(a - 6) ^(wi - Wi) + - Wi) > 0, 

i=l i=l 

or equivalently 

(a - 6)(v(S,) - w(S,)) + 6(v(B,+i) - w(Sf+i)) > 0. 

Since a > b, w{Bg) > v{Bg) and w{Bi^i) > v{Bi^i), the left hand side must be less than or equal to 
zero. We get a contradiction. □ 

Proof of Theorem\16l Suppose without loss of generality that Pi < . . . < Pk and Pi < . . . < Pk- Let r* 
and f* be the max-min fair point in V{g) and 'P{g') respectively. It is clear that g({l}) < • . . < g{{K}), 
5'({1}) < ••■ < 9'{{K}), and both g and g' satisfy the order property. By Theorem [151 we obtain 
r* < . . . < and fl < . .. <f*i^. We want to show 

k k 
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forfc = l,...,/s:-l. 

As in the proof of Prop. [351 ^ bottleneck of the max-min fair vector is of the form {1,2,... 
disregard all constraints except those in the form r(S') < g{S) for 5* = {1, 2, ...,£},£ = 1, . . 
be any integer between 1 and K. By Prop. I36[ the point r* maximizes the region 

ri + ... + rk< 9{{1, 2,...,k}) for fc = 1, . . . , if, 
rk < ru+i for fc 1, . . . , iv: - 1, 

rfc > for fc = 1, . . . , iiT. 

Expressed in terms of matrix, the constraints become r-^A < b^, where 



A = 



111 
1 1 
1 



1 



-1 1 
-1 



and 



= [g({l}),5({l,2}),...,.g({l,2,...,X}),0,...,0] 



'}. We can 
,K. Let k 



Let c denote the column vector [1, . . 
maximal value in the linear program 



1,0,..., 0]"'" in which exactly k components are equal to 1. The 



max rj = ■ r, 

i=i 

r^A < 
r > 

is equal to z^j^i fj- By duality of linear programming, the minimal value of the dual problem, with Sj's 
as the dual variables, 

K 

As > c 
s > 



coincides with J2j=i 

Since P < P', we know that ^({l, . ..,]})> g'({l, . . . , j}) for j = 1, , 
by 5'({1, . . . , we obtain another linear program, 

K 

min^g'({l,...,j})sj. 



,K. If we replace 5({1, ... 



As > c 
s > 0. 



The minimal value equals ■ ^'"^ optimizing a smaller objective function over the same feasible 

region. As a result, we must have X]t=i '^i — ■ ^ 
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